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BACKGROUND OF THE INVENTION 

1 . Field of the Invention 

This invention relates to a method and system for synthesizing a 
circuit representation of a circuit into a new circuit representation having greater 
unateness. 

2. Background Art 

The number of transistors that can be fabricated in a single IC has 
been growing exponentially over the last three decades. A well-known example is 
the Intel series of microprocessors. Intel's first commercial microprocessor, the 
4004, was built with 2,300 transistors in 1971, whereas a recent Intel 
microprocessor, the Pentium III, introduced in 1999 contains 9.5 million transistors. 
The clock frequency of the microprocessors also has dramatically increased from the 
4004's 0. 1 MHZ to the Pentium Ill's 550MHz. The 1998 International Technology 
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Roadmap for Semiconductors developed by the Semiconductor Industry Association 
(SIA) predicts that the transistor count and the clock frequency of ICs will grow 
even faster in the next decade. 

It is thus becoming extremely difficult and time-consuming to design 
5 all the components in a complex IC from scratch, verify their functional and timing 
correctness, and ensure that overall performance requirements are met. To solve 
this challenging problem, a "design reuse" methodology is being widely introduced, 
which integrates large standardized circuit blocks into a single IC called a system 
on a chip (SOC). An SOC integrates a set of predesigned "off-the-shelf" blocks to 
10 build the entire system on a single chip, just as off-the-shelf ICs have been used to 
Jj build a system on a board. The SOC methodology allows IC designers to focus on 

%j the interfaces linking the predesigned blocks. Thus it saves a tremendous amount 

U of time that the designers would have spent creating all the blocks from scratch, and 

H verifying their correctness. For this reason, the SOC design approach is becoming 

15 increasingly popular. 

03 A large portion of the reused blocks in an SOC are intellectual 

jl property (IP) circuits, which are also called cores or virtual components, and are 

*f often provided by third party vendors. The IP providers typically transfer their 

designs to SOC designers in a way that hides the key design details of the IP circuits 
20 and so protect the IP provider's investment in the IP designs. The IP circuits that 

have been developed so far cover numerous functions and support many different 

IC technologies. Additionally, the number and scope of the available IP circuits are 

rapidly growing. 

IP circuits are currently available in three different forms known as 
25 hard, soft, and firm. Hard IP circuits are provided in the form of complete layouts 
that are optimized and verified by the IP providers for a particular IC technology. 
Therefore, hard IP circuits can save time in all SOC design steps, but cannot be re- 
optimized by system designers for other technologies. The intellectual property of 
hard IP circuits includes all the implementation details and is protected by providing 
30 the system designers only with the circuit's high-level behavioral or functional 
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specifications. Soft IP circuits are provided in the form of register-transfer level 
(RTL) descriptions, which define the circuit's behavior using a set of high-level 
blocks. These blocks can be converted by system designers to lower-level designs 
at the gate and physical levels. Thus soft IP circuits can be optimized for a variety 

5 of IC technologies and performance requirements, while they can save SOC design 
time in the high-level design and verification steps. Their RTL designs are 
considered to be the intellectual property contents of soft IP circuits. Finally, firm 
IP circuits are provided as netlists or gate-level descriptions. They allow the system 
designers to optimize the IP circuit's physical design such as cell placement and 

10 routing for various IC technologies. They provide the major advantages of both 
hard and soft IP circuits— they save system design time while allowing the flexibility 
of retargeting the IP circuits at various IC technologies. Both their RTL and gate- 
level designs are considered to be the intellectual property contents of firm IP 
circuits. While hard IP circuits are primarily aimed at ASIC designs, some soft and 

15 firm IP circuits are aimed at SOCs implemented using field programmable gate 
array (FPGA) technology. 

Although the advancement of IC technology allows extremely large 
designs to be integrated in a single chip, today's IC technology presents major 
challenges to the existing design and testing methodologies. For example, testing 
20 requirements for complex digital circuits are becoming increasingly tighter. Using 
traditional processes to synthesize the implementation of digital and other circuits 
often leads to circuits that are either inefficient in meeting testing requirements (i.e., 
unreasonably large test sets, etc.) or cannot satisfy design constraints (i.e., delay, 
area limits, etc.). 

25 The unateness of a circuit has a substantial impact on circuit 

testability and performance. Unate variables of a circuit representation z are 
variables that appear only in complemented or uncomplemented form in z's minimal 
two-level expressions such as sum of products (SOP) or product of sums (POS) 
expressions; binate variables are non-unate. For example, z x = ab + ac + bed is 

30 a minimal SOP expression for the four-variable function z x , in which a and d are 
unate, and b and c are binate. 
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The majority of circuit representations are, in nature, binate, and 
it is often difficult to synthesize these binate functions into a circuit implementation 
that can be efficiently tested for manufacturing defects and operate at very high 
speeds. 

5 For example, a high-speed circuit implementation known as "domino 

logic" requires that the circuit to be implemented be unate. Therefore, binate circuit 
functionality to be implemented in domino logic must be decomposed into unate 
circuit implementations. Similarly, static CMOS logic implementations become 
efficient if unate circuit implementations are extracted from an original binate circuit 
10 prior to implementation. Datapath logic circuits such as adders, subtractors, 
comparators, and ALUs are good examples of applications where carry generation 
functions (i.e., unate functions) are extracted from a larger (often binate) function 
and implemented in high-speed circuit structures. 

Another advantage of unate circuit implementations is their relatively 
15 small universal test set (i.e. , the set of minimal true and maximal false test vectors 
for a function z). These test vectors have the useful property that they can detect 
all multiple stuck-at faults in any implementation of z. Universal test sets guarantee 
a very high coverage of manufacturing defects in a vast range of implementations 
for given circuit functionality. The universal test sets for binate circuit 
20 implementations tend to become excessively large. In addition, the unateness 
property enables the generation of test vectors from the behavioral functions of 
circuits before their implementations are actually executed. 

Existing functional decomposition processes are not suited to the goal 
of decomposing a binate circuit representation into a small set of unate subfunctions. 

25 One existing functional decomposition process called kernel extraction is aimed at 
multi-level logic synthesis. The kernels of a circuit representation/ are defined as 
fs cube-free primary divisors or quotients. For example, the kernels off = (a + 
b + c)(d + e)f + bfg + h include d + e, d + e + g, and a + b + c. This 
decomposition process employs algebraic division operations with the kernels 

30 serving as divisors or quotients to/. Here algebraic division represents /by a logic 
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expression of the form / = pq + r. The kernels of / can be binate, so fs 
subfunctions p, q, and r can also be binate. In addition, kernel extraction often leads 
to an excessive number of subfunctions, and so is not practical for unate 
decomposition. 

5 Another existing decomposition process is Boole-Shannon expansion, 

which represents circuit implementation/by xf x + xf- where f x is the cofactor of 
/ with respect to x. Cofactor f x is defined as the subfunction of / obtained by 
assigning 1 to variable x in/. Boole-Shannon expansion has been widely used for 
binary decision diagram (BDD) construction, technology mapping and field 

10 programmable gate array synthesis. Boole-Shannon expansion is unsuited to the 
goal of obtaining a small set of unate circuit implementations, however, since it may 
only make the child functions f x and f- unate, while always expressing the parent 
function in a binate form. When applied repeatedly, Boole-Shannon expansion can 
also produce an unnecessarily large number of subfunctions, as each is created by 

15 eliminating only one binate variable at a time. 

Finally, a disjoint or disjunctive decomposition represents a boolean 
function /(X, Y) in the form h{g x {X^, g 2 (X 2 ), g^XJ, Y), where X lt X 2 , X k , 
and Y are disjoint (i.e., non-overlapping) variable sets. This decomposition is 
relatively easy to compute, and is sometimes used for logic synthesis problems 
20 where the interconnection cost dominates the circuit's overall cost. It has the 
drawback that many circuit representations cannot be disjointly decomposed, and, 
like Boole-Shannon expansion, it can make the parent function h binate. Thus, the 
disjoint decomposition technique is also not appropriate for our unate decomposition 
goal. 

25 SUMMARY OF THE INVENTION 

One object of the present invention is to provide a method and system 
for efficiently synthesizing a circuit representation into a new circuit representation 
(i.e., circuit implementation) having greater unateness. Notably, those of ordinary 
skill in the relevant art will appreciate that the present invention may be 
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implemented or applied in a variety of circumstances to synthesize circuits beyond 
those discussed, by way of example, in the preceding Background Art. 

One advantage of circuit implementations that possess unateness is 
their relatively small universal test set (i.e. , the set of its minimal true and maximal 
5 false test vectors). Universal test sets guarantee a very high coverage of 
manufacturing defects in a vast range of implementations for a give function. 
Unlike universal test sets for highly unate circuit implementations, the universal test 
sets for largely binate circuits tend to become excessively large. In addition, the 
unateness property enables the generation of test vectors from the behavioral 
10 functions of circuits before their implementations are actually executed. 

Another advantage of unate circuit implementations is their low chip 
area. Yet another advantage of unate circuits is their ability to operate at very high 
speeds. For example, a low-area high-speed circuit implementation known as 
"domino logic" requires that the boolean function to be implemented be unate. 

15 Therefore, binate functions to be implemented in domino logic must be decomposed 
into unate functions prior to implementation. Similarly, static CMOS logic 
implementations become efficient if unate functions are extracted from an original 
binate function prior to implementation. Datapath logic circuits such as adders, 
subtracters, comparators, and ALUs are good examples of applications where carry 

20 generation functions (i.e. , unate functions) are extracted from a larger (often binate) 
function and implemented in high-speed circuit structures. 

To meet these and other objects and advantages of the present 
invention, a method having preferred and alternate embodiments is provided for 
synthesizing a representation of a circuit into a new representation having greater 

25 unateness. The method includes (i) partitioning a circuit representation to obtain 
a representation of at least one sub-circuit, (ii) recursively decomposing the 
representation of the at least one sub-circuit into a sum-of-products or product-of- 
sums representation having greater unateness than the representation of the at least 
one sub-circuit, and (iii) merging the sum-of-products or product-of-sums 

30 representation into the circuit representation to form a new circuit representation. 
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The invented method may additionally include repeating steps (i), (ii) 
and (iii) until a desired level of unateness for the new circuit representation has been 
achieved. 

The invented method may additionally include, for each 
5 decomposition, selecting the sum-of-products or product-of-sums representation 
having fewer binate variables. 

The invented method may additionally include merging common 
expressions of the sum-of-products or product-of-sums representations. 

The invented method may additionally include implementing 
10 algebraic division to merge common unate expressions. 

The invented method may additionally include partitioning the circuit 
representation to obtain a representation of at least one sub-circuit that is highly 
unate. 

The invented method may additionally include implementing a binary 
15 decision diagram to recursively decompose the representation of the at least one 
sub-circuit into the sum-of-products or product-of-sums representation. The binary 
decision diagram may be a zero-suppressed binary decision diagram. 

Additionally, a system having preferred and alternate embodiments 
is provided for synthesizing a circuit representation into a new circuit representation 

20 having greater unateness . The system comprises a computing device configured to 
(i) receive input defining a circuit representation, (ii) partition the circuit 
representation to obtain a representation of at least one sub-circuit, (iii) recursively 
decompose the representation of the at least one sub-circuit into a sum-of-products 
or product-of-sums representation having greater unateness than the representation 

25 of the at least one sub-circuit, (iv) merge the sum-of-products or product-of-sums 
representation into the circuit representation to form the new circuit representation, 
and (v) output the new circuit representation. 
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The computing device may be further configured to receive input 
defining a desired level of unateness for the new circuit representation, and repeat 
steps (ii), (iii) and (iv) until the desired level of unateness is achieved. 

The computing device may be further configured for each 
5 decomposition, select the sum-of-products or product-of-sums representation having 
fewer binate variables. 

The computing device may be further configured to merge common 
expressions of the sum-of-products or product-of-sums representations. 

The computing device may be further configured to implement 
10 algebraic division to merge common expressions. 

The computing device may be additionally configured to partition the 
circuit representation to define a representation of at least one sub-circuit that is 
highly unate. 

The computing device may be additionally configured to implement 
15 a binary decision diagram to recursively decompose the representation of the at least 
one sub-circuit into a sum-of-products or product-of-sums representation. The 
binary decision diagram may be a zero-suppressed binary decision diagram. 

The circuit representation and the new circuit representation may be 
input to the computing device and output from the computing device, respectively, 
20 in a hardware description language such as Verilog or VHDL. 

Additionally, a preferred computer -readable storage embodiment of 
the present invention is provided. In accord with this embodiment, a computer- 
readable storage medium contains computer executable code for instructing one or 
more computers to (i) receive input defining a circuit representation, (ii) partition 
25 the circuit representation to obtain a representation of at least one sub-circuit, (iii) 
recursively decompose the representation of the at least one sub-circuit into a sum- 
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of-products or product-of-sums representation having greater unateness than the 
representation of the at least one sub-circuit, (iv) merge the sum-of-products or 
product-of-sums representation into the circuit representation to form a new circuit 
representation, and (v) output the new circuit representation. 

The computer executable code may additionally instruct the 
computer(s) to receive input defining a desired level of unateness for the new circuit 
representation, and repeat steps (ii), (iii) and (iv) until the desired level of unateness 
is achieved. 

The computer executable code may additionally instruct the 
computer(s) to, for each decomposition, select the sum-of-products or product-of- 
sums representation having fewer binate variables. 

The computer executable code may additionally instruct the 
computer(s) to merge common expressions of the sum-of-products or product-of- 
sums representations. 

The computer executable code may additionally instruct the 
computer (s) to implement algebraic division to merge common expressions. 

The computer executable code may additionally instruct the 
computer(s) to employ a binary decision diagram to recursively decompose the 
representation of the at least one sub-circuit into the sum-of-products or product-of- 
sums representation. 

The above objects and advantages of the present invention are readily 
apparent from the following detailed description of the preferred and alternate 
embodiments, when taken in connection with the accompanying drawings. 
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FIGURES la-lc are functional circuit diagrams illustrating (a) a 
single binate block with three functions, (b) two unate blocks with six functions, and 

(c) three unate blocks with fifteen functions; 

FIGURE 2 is a block representation corresponding to a 
decomposition off(X) = hig^X), g 2 (X), g k (X)) in accordance with the present 
invention; 

FIGURES 3a-3b are (a) a binary tree representing an AND-OR 
decomposition in accordance with the present invention and (b) a corresponding 
block representation in accordance with the present invention; 

FIGURE 4 is a block representation corresponding to a Mevel 
AND/OR tree in accordance with the present invention; 

FIGURE 5 is an example gate-level circuit to be decomposed by the 
computer-implemented process UDSYN in accordance with the present invention; 

FIGURES 6a-6i are (a) a circuit graph G C1 for the circuit of Figure 
5, (b) partition G P1 selected from G C1 , (c) unate decomposition performed on G P1 , 

(d) graph G C12 obtained by merging B 2 and G C1 , (e) partition G P2 selected from G C2 , 
(f) final unate decomposition of G P2 , (g) circuit graph G C3 obtained by merging B 3 
and G C2 , (h) partition G P3 selected from G C3 , and (i) final block representation for 
the circuit of Figure 5 in accordance with the present invention; 

FIGURES 7a-7c illustrate a zero-suppressed binary decision diagram 
(BDD) representing f 4 SOP = ab+ cb : (a) a path representing cube ab; (b) 
a path representing cube cb ; and (c) a zero-suppressed BDD representing 
f/° S - {a + c)(a + b)(b + c) in accordance with the present invention; and 
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FIGURES 8a-8d illustrate four steps in partitioning G C] of Fig. 6a: 
(a) start with a primary input node nlO and select n7; (b) select nl4 and its 
transitive fanin nodes; (c) select nl8 and its transitive fanin nodes; (d) final partition 
G P1 in accordance with the present invention. 



5 DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 



The present invention comprises a method and system having 
preferred and alternate embodiments for efficiently synthesizing a representation of 
a circuit into a new representation having greater unateness. 

For purposes of illustration, the present invention is described in the 
10 context of Bodean function-based digital circuitry. However, application of the 
present invention is not so limited. Notably, the present invention may be applied 
to a variety of circuit representations such as gate-level circuits, PLA 
representation, transistor-level representations, and HDL-based circuits. 

Fig. la defines a circuit representation 10 with three binate outputs 
15 (e.g., x, y and z). Figure lb shows a two-block representation for circuit 10, and 
Fig. lc a three-block representation. The block representations of (b) and (c) are 
obtained by decomposing the functions defining the original single-block circuit 
representation 10. Circuit functionality for each block is represented in a sum-of- 
product (SOP) or product-of-sum (POS) form. Table 1 compares the test 
20 requirements for the three representations of Fig. 1. 
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TABLE 1 



Block 
Representation 


Test Requirements 


Implementation Flexibility 


No. of Binate 
Variables for 
all Functions 


Universal 
Test Set 
Size 


No. of 
Block 
Functions 


Area of an Example 
Synthesized 
Circuit 


Figure la 


7 


64 


3 


28 


Figure lb 


0 


33 


6 


29 


Figure lc 


0 


47 


12 


50 



In the case of Fig. la, function z is binate for all six variables, so its 
universal test set consists of all 64 possible test vectors. The decomposed designs 
of Figs, lb and c require smaller universal test sets (33 and 47, respectively), since 
all their internal blocks are unate. Although the difference in test set size is minor 
in this small example, it tends to be significant for larger circuits. 

Circuits that do not have natural block representations are often 
implemented by logic synthesis systems, while those with natural block 
representations are often implemented manually. Although the present invention is 
not limited to the former case, we assume that the target circuit is of the former type 
in order to best describe the present invention. 

For a given circuit, the block representations created in accordance 
with the present invention can be considered as design constraints. In other words, 
the boundaries of the blocks serve as a high-level structural design constraint that 
must be satisfied by low-level implementations of the circuit such as gate-level or 
transistor-level implementations. Notably, these design constraints tend to restrict 
implementation flexibility. 

The outputs of any blocks in a circuit C's block representation S B 
define the block functions. Circuit C's implementation flexibility can be roughly 
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measured by the number of block functions that C employs. This follows from the 
fact that a large number of block functions in S B implies that the functions 
themselves are small. Block functions in different blocks cannot be merged, so 
implementations of small block functions are generally less flexible than large ones. 
Thus the fewer the block functions in S B , the higher the implementation flexibility 
of C. 

For example, Table 1 compares the implementation flexibility of the 
block representations in Fig. 1. Figure lb has three more block functions than Fig. 
la, while Fig. lchas 12 more. Whereas Fig. la has full implementation flexibility, 
the other block representations have limited flexibility. The last column of Table 1 
compares the area of some example implementations synthesized by a commercial 
synthesis system Synopsys® Design Compiler with the goal of reducing area. The 
area is calculated from the relative gate areas defined in the Synopsys cell library; 
for example, inverter = 1, AND2 = 2, AND4 = 3, OR2 = 2, and OR4 = 3. In 
summary, this example suggests that lower implementation flexibility often leads 
to poor implementations in terms of circuit area. 

To permit a broad range of implementation styles, the invented 
synthesis process attempts to decompose a binate function into as small a set of 
unate subfunctions as possible. In general, a decomposition of function /can be 
expressed as: 

f(X) = h( gl (X),g 2 (X), ..., gk (X)) 

Let/be the root function, subfunction h the parent function, and each 
subfunction g, a child function. X = {x u x 2 , x n } denotes the support set off and 
each g;. A decomposition is called a unate decomposition if all the subfunctions h 
and g Lk inf(X) are unate. A decomposition of the form/(X) transforms a single- 
block model of function/into a two-block model with h defining one block and gl . k 
the other; see Figures 2 and 3. 
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In accordance with the present invention, a preferred method for 
synthesizing circuits utilizing unateness properties of boolean functions involves 
recursive OR and AND decompositions of a circuit output function / and its 
subfunctions. The OR and AND decompositions represent f(X) by h(g^(X), g 2 (X)), 
5 where h has the form g Y + g 2 and g^ g 2 , respectively. We obtain an OR 
decomposition off from an SOP form off, and an AND decomposition from a POS 
form. 



For example, consider a binate function f whose SOP form is 
ab + ac + ca . A possible OR decomposition of f is h = g t + g 2 with g x = ac 
10 and g 2 = ab + ca . This decomposition makes all the subfunctions h, g u and g 2 
unate, and so is a unate decomposition. Now consider /j's POS form 
(a + c)(a + b + c) . We can obtain directly from this form an AND 
decomposition with subfunctions h = g 2 , g x = a + c, and g 2 = a + b + c . 
This is also a unate decomposition. 



15 A single OR or AND decomposition of a large binate function may 

not lead to a unate decomposition. However, a sequence of OR or AND 
decompositions applied to /recursively always produces a unate decomposition for 
any function /. The general form of such a sequence with k levels is 



f=hKg^g 2 J ) 

20 ft 1 = h 2 ( §1 \ g 2 2 ) 

£; k l = h\g x k , g 2 *) 



where h! and g{ denote a parent function and a sub function, respectively, produced 
by the j'-th level decomposition. Parent function h> can be either AND or OR. A 
25 Ar-level sequence of AND and OR decompositions forms a binary AND-OR tree, 
where the internal nodes represent AND or OR operations, while the leaf nodes 
denote unate subfunctions that need no further decomposition. 
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An arbitrary sequence of AND and OR decompositions can lead to 
an excessive number of sub functions. To reduce this number, we restrict our 
attention to sequences of the following type, which we refer to as unate AND-OR 
decompositions . 

5 f=h\g u \ 8b 1 ) 

g„ l = h 2 (g u \ g b 2 ) 

8b' 1 = hk (8u, g b k ) 

As in the general case, h> is either AND or OR, but the final g b k and 
every gj are unate, while every g b j except the final one is either unate or binate. 
This decomposition can also be represented in the compact factored form 

f=h\g u \ h\g 2 , h\gi, ... , fr l (gf\tt(g a k ,g b *))...))) (1) 
as well as the general form 

15 m = Kg u \g u \ ...,g u k , 8b k ). (2) 

Comparing (1) with (2), we see that the parent function h in (2) is 
composed of the AND and OR subfunctions h\ h 2 , ...,h k only, and so is always 
unate. 

Figure 3a shows a binary AND-OR tree corresponding to the unate 

20 decomposition 

/= (gu-igu 1 * (Su* + (gS-(g u 5 ■&?))))) 

obtained by a unate AND-OR decomposition with 5 levels. The internal nodes 30a- 
30e in Fig. 3a represent AND or OR operations, while the leaf nodes 32a-32f at the 
bottom represent unate subfunctions. Figure 3b depicts the block representation 
25 corresponding to Figure 3a. Bj defines an AND-OR tree network that implements 
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the function h. B 2 is a network of undefined internal representation that implements 
the unate subfunctions g u \ g u 2 , g u \ g u \ g u \ and g b 5 . 

In general, we obtain the foregoing kind of unate AND-OR 
decomposition for/as follows: first decompose/ into g u l and g b l using an AND or 
5 OR operation that makes g u l unate; then repeatedly decompose gj into gj +1 and g b J+1 
in a similar way, until g b J+1 becomes unate. This must eventually happen, because 
a gj of a single product or sum term is unate. In practice, the AND-OR 
decomposition process often terminates with a final g b J consisting of a relatively 
large unate function. 

CI 10 As noted above, the global parent function h(g u \ g u k , g b ) in (2) 

JI is unate. Thus, the final result of a £-level AND-OR decomposition is a set of k + 

W 2 unate subfunctions g u l :k , g b k , and h{g u l , . . . , g u \ g b k ) . 

f: Notably, an important goal of the block synthesis method shown in 

Figure 3 is to find an AND-OR decomposition of a given function /using as few 
15 subfunctions as possible. In addition, it is preferred that each gj be selected in a 
=3" manner that makes the resulting gj highly unate. This selection often leads to a 

Pi unate decomposition involving few subfunctions. Also such a gj can be relatively 

^ easily derived from a standard SOP or POS form. 

Each level of a unate AND-OR decomposition is defined by either 
20 an AND or OR operation. How we select the operation at each level has a large 
impact on the final result, as we show with the following example. 

Consider^ = a © b © c whose SOP and POS forms are given below. 



f 2 = abc+ abc+ abc+ abc (3) 
f 2 POS = (a+ b+ c)(a+ b+ c)(a+ b+ c)(a+ b+ c) (4) 
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OR decompositions are derived from (3), and AND decompositions 



are derived from (4). Suppose we select an OR operation in every level of the 
decomposition. A possible result is: 



Note that in this particular example, the unate decomposition is 
completed when the final g b k is a product term, and the resulting g 1 subfunctions 
correspond to each of the product terms in (3). 

Next, suppose we select an AND operation in every level. A 
15 possible result is the unate decomposition^ = g u l -(g u 2 (g„ 3 -g b 3 )), which involves 
five subfunctions: 



fi = gu + igu + (gu* + g b 3 )) 



5 



which involves five unate subfunctions 



gu = ab ° 
gu = ab ° 



10 



g b 2 = abc ; and 



g> 



u 



a\ b+ c ; 



2 



a+ b+ c ; 



g, 



u 




a+ b+ c ; 



20 



g b 3 = a + b+ c ; and 



h(gu 'gu >gu 'gb ) • 
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Notably, a unate decomposition of f 2 can be obtained involving fewer 
subfunctions if AND and OR operations are mixed as follows. Suppose we select 
an OR operation in the first level and an AND operation in the second level. The 
OR operation decomposes (3) into f 2 l = g u ' 4- g b \ where g J = ale and 
5 = abc+abc+ abc . To apply an AND operation to g^, 1 , we use g b l 's POS form 

(a + b)(a + c)(b + c)(a + b + c) . Then an AND operation leads to g b = g u z - g b 2 , where 
for example, g u 2 = (a+b+c) and g b 2 = (a + V){a + c)(b + c) . Since g b z is unate, 
the unate decomposition is complete. Note that unlike the previous cases, the final 
gb here contains more than one term. The third unate decomposition of/ 2 is 

10 

I fi = Su + (gu ■ g b Z ) = abc + (a + b + c) ■ (a + b)(a + c)(b + c) 



which involves only four subfunctions, one less than the first and second cases, 
where we selected three AND and three OR operations, respectively. This example 
shows that how we select the AND-OR operation in each level of the AND-OR 
15 decomposition is very important. 

Often, there are many possible AND and OR decompositions in each 
level. This implies the existence of a large number of possible unate AND-OR 
decompositions. For example, if an SOP form of gj contains m product terms, we 
can partition these terms into two groups defining gj +1 and g b J+1 in 2 m different 

20 ways. At each level, either an AND or OR operation can be chosen, so the number 
of possible &-level unate AND-OR decompositions is 2 m+k . Thus, finding a unate 
AND-OR decomposition of a large function / involving a minimal number of 
subfunctions is often impractical. We therefore introduce Unate-Decomp, a 
heuristic process that systematically selects AND or OR decompositions at each 

25 recursion level, and produces a final unate AND-OR decomposition containing 
relatively few subfunctions. 

Unate-Decomp represents a function /and all its subfunctions in both 
SOP and POS forms. To produce an AND (OR) decomposition off, it selects a set 
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S of product terms (sum terms) from the SOP (POS) form of /, so that 5 constitutes 
a unate subfunction g u . The rest of the product terms (sum terms) of / define 
subfunction g b . Unate-Decomp then represents g b by both SOP and POS forms, 
which it uses to produce an OR and AND decomposition at the next recursion level. 
5 To represent SOP and POS forms efficiently, binary decision diagrams can be 
employed. 

To decompose / into as few unate subfunctions as possible, Unate- 
Decomp produces each gj in a way that reduces the number of binate variables in 
g b 3 . In the case of multiple- function circuits, Unate-Decomp first decomposes each 
10 output function f t using the method described above. It then merges common 
subfunctions of different functions f t and fj to reduce the total number of 
subfunctions . 



Notably, representing large circuits directly by two-level expressions 
is often inefficient. To handle such cases efficiently, a preferred process first 

15 partitions a given circuit, and then performs decomposition on each partition. For 
example, one partition is created for each process of Unate-Decomp parent function 
h. The resulting decomposition is then merged into the rest of the circuit. Then the 
next partition is created from the merged circuit, and the next process of Unate- 
Decomp is conducted. This process is repeated until no more partitioning is 

20 necessary. 



Preferably, each decomposition step and circuit partition are selected 
in a way that produces a small number of highly unate subfunctions. Consequently, 
the resulting block representations tend to have a high level of implementation 
flexibility. 

25 To decompose /into as few unate subfunctions as possible, Unate- 

Decomp produces each gj in a way that reduces the number of binate variables in 
g b j . For example, consider the following function/: 
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f 3 = abc + ad + ae + bcf + bdf + be + cdf + ce + abdf 
+ abe + acdf + cdf + cd + de 

Suppose we decompose / 3 into g u + g b . Table 2 shows some 
possible ways of doing this and the number of binate variables in the resulting g b . 



TABLE 2 



gu 




No. of 
Binate 
Variables in 

g b 


cd + acd f + bd f 


ad+abe+bcf+cd f + abc+ abd f 
+ ae+ de + cd f + be+ ce 


6 


abd f + ce 


ad+abe+bcf + cd f + abc+ ae 
+de+cdf+be+cd+acdf+bdf 


5 


bcf + bdf + cdf + be + 
ce + abc + ae 


cd + acdf + abe + cdf 
+ ad + abd f + de 


3 


abc + ad + ae + bcf + b 
df + be + cdf + ce 


abdf + abe+ acdf 
+ cdf + cd + de 


2 



While the first OR decomposition produces g b with six binate 
variables, the last OR decomposition produces g b with only two binate variables, 
and so is selected. 
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In each level of the decomposition process, Unate-Decomp produces 
a pair of AND and OR decompositions using a special form of cofactor operation 
called Subset. The Subset operation for a literal extracts a subset S of product 
(sum) terms of a given SOP (POS) form by eliminating terms that contain l t . For 
5 example, applying Subset to the SOP form 

abc + acd + abd 

for literal 

a 

yields 

10 S - acd + abd 

Unate-Decomp systematically computes S for a set of binate literals 
{1} so that 5 is unate and the set of other terms is highly unate. Then S defines gj, 
and the other product terms define gj. 

After a unate decomposition is formed, Unate-Decomp constructs two 
15 blocks from an AND-OR tree representing the decomposition; see Fig. 4. To 
ensure that all the block functions are unate, we place in block B x all the nodes 
representing the subfunctions g u h k and g b k , which correspond to the leaf nodes in the 
AND-OR decomposition tree. We place in block B 2 all the other nodes, which 
represent AND and OR operations and together form the function h. 

20 In the preceding description, we focused on decomposing a single 

function. In the case of multiple-function circuits, Unate-Decomp first decomposes 
each output function^ using the method described above. It then merges common 
subfunctions of different functions / t and to reduce the total number of 
subfunctions. Algebraic division operations are often employed by logic synthesis 

25 techniques to efficiently combine common expressions. These operations can be 
easily applied to the results of our AND-OR decompositions, and often reduce the 
number of subfunctions significantly. Notably, Unate-Decomp incorporates 
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algebraic division in such a manner that two different functions share the divisor of 
each division. 

Based on the unate decomposition concept described above, we 
introduce a computer-implemented synthesis program in accordance with the present 
5 invention called Unate Decomposition Synthesis (UDSYN). Representing large 
circuits directly by two-level expressions is often inefficient. To handle such cases 
efficiently, UDSYN first partitions the given circuit, and then performs 
decomposition on each partition, as generally described above. 

10 Table 3 contains one embodiment of a pseudocode representation of 

UDSYN in accordance with the present invention. Notably, it is understood by those 
of ordinary skill in the art that different computer programs and program 
arrangements can be implemented to support and execute the overall function of 
UDSYN. 

15 TABLE 3 



Embodiment of process UDSYN (Verilog-input) 




1 : G c : = Build-Circuit-Graph(y erilog- input) ; 




2: while (G c * 0) begin 




3: G P := UD-Partition(G c ); 


/* G P is a graph representing a partitioned block */ 


4: G c := G c - G P ; 


/* Remove nodes in G P from G c *l 


5: for each output node n R in G P begin 




6: Build-ZSBDD (n R , G P ); 


1* Create SOP and its complement for n R *l 


7: end; 




8: (5„ B l+i ) := Unate-Decomp(Gp); 


/* B : and correspond to B x and B 2 of Fig. 7 */ 


9: G c := G c uB, +l ; i : = i + 1; 


/* Insert notes in S, +1 into G c */ 


10: end; 




11: Verilog-output := Interconnect-Blocks{{B,}); 


/* Verilog-output is the final block representation */ 


12: return Verilog-output; 





UDSYN takes an input circuit in a form such as a Verilog 
30 specification whose internal elements can be either functional descriptions or gates. 
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First, UDSYN builds a circuit graph G c whose nodes represent the internal elements. 
It then creates a partition G P of G c using UD-Partition(G c ) , and removes nodes in 
G P from G c . The output functions of G P are represented in SOP and POS forms. 
The process Unate-Decomp(G P ) performs unate AND-OR decompositions on the 
5 output nodes in G P , and constructs decomposed blocks B : and B 2 as in Fig. 4. 
Blocks B 1 and B 2 created from the i-\h partition G Pi are denoted by B { and B i+1 , 
respectively. Step 9 modifies G c by inserting all nodes of B 2 into G c . UDSYN 
repeats the above steps until all nodes in G c are removed. It then constructs a 
hardware description language (e.g. Verilog, VHDL, etc.) output file by specifying 
10 the interconnections among all blocks B { . 

We illustrate UDSYN using a gate-level circuit of Fig. 5 as input. 
Figures 6a to i show intermediate results of steps 2 to 10 in Table 3. Figure 6a 
shows the circuit graph G C1 for the circuit of Fig. 5; each node in G C1 corresponds 
to a gate in the circuit. UD-Partition(G c ) creates a partition G P1 starting from the 

15 primary inputs of G C1 . The shading in G C1 indicates nodes that are selected by UD- 
Partition(G c ) , and constitute G P1 . Figure 6b represents G P1 by a rectangle. All 
nodes in G P1 are removed from G C1 , and are merged into SOP and POS forms by 
steps 5 to 8. These SOP and POS forms are decomposed by step 8 into unate 
sub functions; these sub functions are grouped into two blocks B x and B 2 as in Fig. 

20 4. As Fig. 6c shows, B r consists of seven subfunctions and B 2 consists of three 
subfunctions. We create a new circuit graph G C2 by merging B 2 and G C] as shown 
in Fig. 6d. Returning to step 3, UD-Partition(G c ) selects some nodes (shaded) in 
G C2 and creates a new partition G n , which is represented by a rectangle in Fig. 6e. 
Then step 8 decomposes G^ into blocks B 2 and B 3 appearing in Fig. 6f. Step 9 

25 merges B 3 into a new circuit graph G C3 as in the preceding steps; see Fig. 6g. 
Figure 6h shows a new partition G P3 constructed from G C3 . By repeating this 
process, we finally obtain the decomposed block representation of Fig. 6i consisting 
of five blocks B 15 . The output functions of these blocks are described by Verilog 
code in equation form. If UD-Partition(G c ) constructs k partitions, UDSYN 

30 produces a total of k + 1 blocks. 
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Step 6 of Table 3 uses a type of binary decision diagram (BDD) 
called a zero-suppressed BDD (ZSBDD) to represent the SOP and POS forms of 
functions. Although other forms of BDD can be used, we limit our attention in this 
description to ZSBDDs for the sake of presentation. A ZSBDD of a function / is 
a graph whose paths denote the product terms (cubes) in an SOP form of /. 
UDSYN uses two ZSBDDs to represent a pair of SOP and POS forms for the 
internal function of each node in an AND-OR tree like that in Fig. 4. Thus an 
AND-OR tree with n nodes is represented by In individual ZSBDDs, each of which 
is linked to the corresponding node in the tree. 

For example, Figs, la and b show a ZSBDD representing 
f i sop = ab+ct . The internal nodes (circles) in Figs, la and b denote the literals 
appearing mf 4 SOP . The terminal or leaf nodes (rectangles) denote the output values 
that f 4 SOP generates when its literals are set to the values specified on the edges. The 
name "zero-suppressed" stems from the property that all nodes in a ZSBDD whose 
1-edge points to the 1 -terminal node are eliminated. Every path from the root to the 
1 -terminal node represents a cube inf/ op . For example, the path highlighted by the 
dashed line in Fig. la represents cube ab, while the one highlighted in Fig. lb 
represents cube cb . Although ZSBDDs can represent only SOP forms directly, 
POS forms can also be handled by their complemented form. 

For example, consider the POS expression 

f/ os = (a+ c)(a + b)(b + c) 

having the following complement 

f A pos = ac+ab+bc 

Figure 7c shows a ZSBDD that represents f/ os , where every path 
from the root to the 1 -terminal node represents a sum term mf/ os with their literals 
complemented. In this way, we can represent both SOP and POS forms using 
ZSBDDs. 
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ZSBDDs have been shown to represent large functions efficiently. 
This is due to the fact that small ZSBDDs can contain a large number of paths, so 
a large number of cubes can be often represented by a compact ZSBDD. ZSBDDs 
also support fast manipulation of sets of cubes such as finding subsets, computing 
5 the intersection of cube sets, and performing a type of division operation on cube 
sets. Utilizing these features, we implement unate AND-OR decomposition and 
division processes that act directly on ZSBDDs. 

UD-Partition(G c ) creates a partition of the input circuit in a way that 
makes the functions defining the partition highly unate, while meeting a specified 
10 partition size limit. Partitions created in this way simplify the unate decomposition 
process. One pseudo-code embodiment of UD-Partition appears in Table 4. 
Notably, it is understood by those of ordinary skill in the art that a variety of 
different computer programs and program arrangements can be implemented to 
support and execute the inventive function of UD-Partition. 



15 



TABLE 4 



20 



25 



30 



Embodiment of process i7D-Partition(G c ) 
1 : for each n c in G c in level order begin 
2: for each fan in node n, of n c begin 
3: if (n, is a primary input of G c ) then 

4: S supp {n c ) : = S supp (n c ) u n,; 

5: else 

7: Calculate-Path-Count{S supp {n c )); 
8: end; 

k 

9: N BV (n c ) = ^ Binateis, , n c ) 

10: end; 

11: while (G c *0) begin 

12: n m := Select-Node-of-Min-N BV (G c ); 

13: S N : = nodes in n,„'s fan-in cone in G c ; 

14: if (I/0-count(G P u S N ) < threshold) then 

1 5 : Gp '. = Gp u S N ', Gq '. = G c — G P \ 

16: else break; 

17: end; 

18: return (G c , G P ); 



/* S SU p P ( n c) is n c' s support set */ 



/* Add path count of fan-in nodes to n c *l 



/* Binate(s„ n c ) = 1 if s, is binate for n c */ 



/* n m is to be included in the partition */ 

/* Gp is the graph for the partition */ 

/* Add «„, and its transitive fan-in nodes to G P *l 

I* Discard the candidate node n m */ 
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Steps 1 to 10 compute the number of binate support variables of each 
node in the circuit graph G c . Steps 11 to 17 create the current partition G P by 
selecting nodes in G c that have a minimal number of binate variables. 

A node n c in G c is unate with respect to a primary input s u if all paths 
between n c and s { have the same inversion parity; otherwise, n c is binate with respect 
to s,. To determine the inversion parity of the paths, we calculate the number of 
paths from the primary inputs to each node n c in G c . Let p^is,, n c ) and p^fa, n c ) 
be the number of paths from a primary input s, to a node n c whose inversion parity 
is even and odd, respectively. Steps 3 to 7 find the set S supp (n c ) of support variables 
for each node n c . For n c and its fanin nodes n„ Calculate-Path-Count obtains p^s^ 
n c ) and p^s,, n c ) by recursively computing 

Aven^i, «c) = /Wnfo, « c ) + ^evenC?,, «i) 

PodM, « c ) = Poddi^, « c ) + AxidOi, 
if the inversion parity from n K to n c is even; otherwise, it computes 

/'even 

(5;, Tl c ) = /7 e ven(X, «c) + AxidOi, «,) 
PoddOi, « c ) = Poddfap K C ) + AvenOi, ",)• 

The binary function Binate(s { , n c ) produces 1 (0), if a node n c is 
binate (unate) with respect to its support variable s u and is computed by 

Binate{s { , n c ) = 0, if p^fr, n c ) = 0 or Poii { Sl , n c ) = 0 
Binate(s { , « c ) = 1, otherwise 

The number N m (n c ) of binate variables of node n c with £ variables is 

defined as 
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The intuition behind using N B ^n^) to guide the partitioning stems 
from the fact that the more binate the node « c , the more difficult the decomposition 
process for n c tends to be. Steps 1 to 10 traverse every node only once, which has 
complexity 0(N). They propagate p even (s h « c ) and p^fa, n c ) for every s { for a n c to 
5 the nodes in the transitive fanout of n c , which also accounts for complexity 0(N). 
Hence the overall complexity of computing N BV (n c ) for all nodes in G c is 0(1^). 

For example, Table 5 shows the calculation of N m (n,.) for every node 

in Fig. 6a. 



TABLE 5 



; 1 O 


XT ^ J „ 

Node 


No. of paths from n c s support variable j, 


No. of binate 
variables 




n5 


(a, 0, 1), (6, 0, 1) 


0 




nlO 


(ft, 0, 1), (d, 0, 1) 


0 




n3 


(f, 1, 0), {d, 1, 0) 


0 


tti 15 


n\2 


(g, 0, 1) 


0 




nl3 


(A, 1, 0), (j\ 1, 0) 


0 


3 '3 


n8 


(c, 0, 1), (ft, 1, 0), (rf, 1, 0) 


0 




nl 


(£, 1, 0), if, 0, 1), <rf, 0, 1) 


0 




nl4 


(h, 1, 0), (i, 1, 0), (f, 1, 0), (rf, 1, 0) 


0 


20 


n9 


(c, 1, 0), (ft, 0, 1), (d, 0, 1) 


0 




nil 


(ft, 1, 1), id, 1, 1), (c, 1, 0), (e, 0, 1) 


2 




n4 


(c, 1, 1), (ft, 1, 1), (d, 2,2), (g, 1, 1), <jf, 1,1) 


5 




ril 


(c, 0, 1), (ft, 2, 0), (d, 1, 0), (a, 1, 0) 


0 




n6 


(a, 1, 1), (ft, 2, 2), (c, 1, 1), (d, 1, 1) 


4 


25 


«18 


(c, 1, 1), (ft, 1, 1), (</, 3, 2), (g, 1, 1), CT, 2, 1), (A, 1, 0), (/, 1, 
0) 


5 




«15 


(c, 2, 2), (ft, 4, 4), (d, 3, 3), (a, 1, 1), (e, 1, 1) 


5 




nl 


(f, 0, 1), (d, 1, 2), (a, 1, 1), (ft, 2, 2), (c, 1, 1) 


4 
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Node 


No. of paths from n's support variable s, 

X L X X I 


No. of binate 
variables 

N BV (n c ) 


nl6 


(f, 3, 3), (d, 1, 7), (a, 2, 2), (b, 6, 6), (c, 4, 4), (g, 2, 2) 


6 


nil 


(g, 5, 5), (f, 7, 7), (d, 15, 15), (a, 4, 4), (&, 12, 12), (c, 8, 8) 


6 



The second column lists p^Si, n c ) and p^is,, n c ) computed for each 
node n c and all its support variables. The last column gives N BV (n c ). For example, 
for n c = nil, Binate(b, nil) = 1, Binate(d, nil) = 1, Binate(c, nil) = 0, and 
Binate(e, nil) = 0. Thus N BV (n c ) = 1 + 1+ 0 + 0 = 2. 

After N BV (n c ) is computed for every n c in G c , UD-Partition selects 
from G c a node n m of minimal N BW (n c ) starting from a primary input of G c . It then 
inserts into G P all non-partitioned nodes in the transitive fanin region of n m . This 
process is repeated until the size of G P exceeds a threshold equal to the maximum 
number of G P 's I/O lines. By limiting the partition size in this way, we can prevent 
ZSBDDs from exploding for large circuits, while producing a partition with highly 
unate output functions. 

Figure 8 illustrates how we partition G c of Fig. 6a. Suppose we limit 
Gp's I/O lines to seven inputs and six outputs, that is, we set the threshold to 7/6. 
The N BV (n c ) values calculated in Table 5 are shown next to each node n c in Fig. 8. 
Figures 8a to d indicate the current G P created in each iteration by shading, and 
newly selected nodes by thick circles. The first n m is selected from the candidate 
nodes n3, n5, n8, nlO, nil, nl2, and nl3, which are adjacent to the primary inputs. 
We select «3 whose N BV (n c ) has the minimum value 0, and add it to G P ; see Fig. 8a. 
The next search begins from n3 and selects n3's fanout node nl whose N BW (n c ) = 
0. We then select all nodes in the transitive fanin region of nl; Fig. 8a indicates 
these selected nodes by a dashed line. Figure 8b shows the current G P consisting 
of rii, nil, and nl. We then select nl4 over nl, nA, and nl7, and then select nl4's 
fanin node nl3; the newly selected nodes are again indicated by a dashed line in 
Fig. 8b. We next select nil over nl8, nl, nA, but nil leads to a partition with 
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seven outputs, one greater than the limit six. Hence we select nl8 instead which has 
the next smallest N BV (n c ). We then select nodes in «18's transitive fanin region; see 
Fig. 8c. At this point, the number of I/O lines of G P equals the threshold 7/6, so 
the partitioning is done. Figure 8 d indicates the final G P by shading. 



10 



Since UD-Partition selects nodes with fewer binate variables first, it 
often leads to a partition where many output functions are already unate and so 
require no further decomposition. For example, in G P of Fig. 8d, four output 
functions at n3, ril, n8, and n 10 are unate. Figure 6c shows a unate decomposition 
of this G P , where nodes g3, gl, g8, and #10 in D x correspond to these four unate 
functions, and so are not decomposed. 



Next we describe Unate-Decomp(G) which systematically applies 
unate AND-OR decomposition operations to a circuit partition. See Table 6 for one 
pseudo-code embodiment of Unate-Decomp(G) . 



TABLE 6 



Embodiment of process Unate-Decomp(G) 




1: 


S B : = G's binate function nodes; /* S B 


stores nodes of binate functions to be decomposed */ 


2: 


while (5 B * 0) begin 




3: 


for each node n B in S B begin 




4: 


(n u , n b ) : = ANDOR-OneLevel(G, n B ); 


/*«„ (n b ) points to subfunction g u (g b ) in (3.4) */ 


5: 


S D := S D u {n u , n b ); 


/* S D stores candidate divisor nodes */ 


6: 


end; 




7: 


for each node n d in 5^ begin 




8: 


for each node ^in S B - S D begin 


/* n f is a candidate dividend node */ 


9: 


(n q , n r ) := Division(n f , n d ); 


/* n q is the quotient and n r is the remainder */ 


10 


: if {N^tif) < N BV Qi 9 ) + N BV (n r )) then 




11 


: Reverse the division; 




12 


end; 




13 


: end; 




14 


: S B := 0; S D := 0; 




15 


for each node n t in G begin 


/* Find new nodes to be decomposed */ 


16 


: if (N BY (n,) > threshold) 




17 


: S B : = S B u »,; 


/* «, exceeds the threshold */ 


18 


end; 




19 


: end; 




20 


: return G; 





15 



20 



25 



30 



35 
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Graph G initially contains the nodes in the current partition. Steps 
3 to 6 perform a level of AND-OR decomposition on every binate node n B in G. 
Then, steps 7 to 13 perform division operations on every binate node in G by 
treating as divisors child nodes created by the AND-OR decompositions. N BW (n-) 
denotes the number of binate variables in the subfunction at node n x in G. If a 
division reduces N BV (n-X it is accepted; otherwise, it is discarded. The above 
process is repeated until all nodes in G become unate. For some large binate 
functions, forcing all nodes to be unate leads to an excessive number of 
subfunctions. We therefore stop decomposing a node n { if A^v^i) becomes less than 
a certain threshold. This threshold is chosen to yield a small set of subfunctions at 
the cost of lower unateness. Thus the threshold allows us to trade the level of 
unateness for a higher implementation flexibility of the block representation. 

Table 7 contains a pseudo-code embodiment of the computer- 
implemented process ANDOR-OneLevel{G, n B ), which implements one level of the 
AND-OR decomposition technique described earlier. 

TABLE 7 

Embodiment of process ANDOR-OneLevel(G, n B ) 

SOP SOP 

1: (Su >gb )■= Find-Unate-Cube-Set (SOP(n B )); /* OR decomposition */ 

CPOS CPOS 

2: (g u >b ) : = Find-Unate-Cube-Set (Inv(SOP(n B ))); /* AND decomposition */ 

SOP CPOS 

3: itiN By (g b )< N BV (g b )) then 

4: Replace-Node (n B , NewNodes (h sop ,gJ° F ,g b S °" )) ; /* Replace n B in G by the new nodes */ 
5 : return (Node(g u *°* ),Node(g b )) ; 

6: else 

7: Replace-Node (n B , NewNodes (h pos ,Inv{g"° s ),/«v(g i <= ™ ))) ; /* ln\(g) complements g */ 
8: return (Nodeilnvig™ )\Node{Inv(g b """))) ; 



The process Find-Unate-Cube-Set(SOP(n B )) finds a set of unate cubes 
(product terms) from an SOP representation SOP(n B ) for node n B . This operation 
forms an OR decomposition. An AND decomposition is obtained by 
complementing the input SOP(n B ) and the outputs of Find-Unate-Cube-Set, 
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respectively. This enables ZSBDDs to handle both AND and OR decompositions, 
although ZSBDDs can only represent SOP forms directly. 

Table 8 contains a pseudo-code embodiment of the computer- 
implemented process Find-Unate-Cube-Set(ISOP) . 



TABLE 8 





Embodiment of process Find-Unate-Cube-Set (ISOP) /* ISOP is the initial SOP form */ 




1: 


S besl : = SOP : = ISOP; 




2: 


while (N BY (SOP)> threshold) begin 




3: 


for each literal /, for binate variables in SOP repeat 


4: 


S, : = Subset(SOP, I,); /* Remove all cubes containing /, */ 




5: 


if (N^+N^ISOP - S,)<N BV (S bei ,)+N BV (ISOP - S best )) then 




6: 


Sbest - = S t ; 




7: 


end; 




8: 


SOP:= S besl ; 


h 15 


9: 


end; 




10: 






11: 


&:=/50P-5 to ; 


3 2 


12: 


return (g u , g b ); 



Find-Unate-Cube-Set{ISOP) derives a cube set S from/s initial SOP 
20 form /SOP so that S meets the threshold on N BV (S). As a result, S defines unate 
subtraction g/ in (1), while ISOP - S defines g b k . As discussed earlier, an exact 
method to find an optimum AND-OR decomposition of an w-term ISOP must 
examine up to 2 m possible AND decompositions. To avoid this and derive S 
efficiently, a type of cofactor operation is implemented which can simultaneously 
25 extract from ISOP multiple cubes (product terms) with a common property. This 
operation, denoted by Subset{SOP, removes from SOP all cubes that contain 
literal / ; . Thus /; does not appear in the resulting SOP form S„ while I may appear 
in S,; hence 5, is unate with respect to / ; . 

For example, consider a function f 5 whose SOP form is 
30 SOP = abc + acd + ad + bda + ed 
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Subset(SOP, d) removes cubes ad and bda which contain literal 

d , and yield 

S, = abc + acd + ed 

The basic concept is to apply Subset(SOP, lj to a set of binate literals 
5 in a way that makes both S, and SOP - S { highly unate. We found that for a binate 
k, Subset(SOP, /;) often eliminates from SOP not only /, but also other binate 
literals. Hence we can often obtain a highly unate cube set S by repeating only a 
few steps of Subset{SOP, /,). The inner loop (steps 3 to 8) of Table 8 performs 
Subset{SOP, I) for all binate literals /„ and selects S besl , (i.e., the 5, having the 
10 minimum A' BV (5 I ) + N BV (ISOP - S-J). The outer loop (steps 2 to 9) repeats this 
process recursively with 5 best in place of SOP until A^vC^best) becomes less than 
threshold. The final S hest defines g u , while ISOP-S best defines g b . 

To illustrate, consider an initial SOP form 

ISOP = abc + acd + ad+ bda + ec + abde + ed 

15 Suppose that the threshold of A^ BV is 0. Table 9 shows each step of the outer loop 
in its first iteration with ISOP assigned to SOP. 



TABLE 9 



20 



SOP = ISOP = abc + acd + ad + bda + ec + abde + ed 


Binate 
literal l t 


Si 


ISOP - s t 


N BV (S,) + 
N^ISOP-S,) 


a 


be + abde + ed 


abc + acd + ad + bda 


3 


a 


abc + acd + ad + bda 
+ bc+ ed 


abde 


3 


b 


ad + bda + ed+bc 


abc + acd + abde 


2 


b 


abc + acd + ad + ed + abde 


bda+ be 


2 


c 


abc + acd + ad+ bda + abde 


be 


2 
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SOP = ISOP = abc + acd + ad + bda + ec + abde + ed 


Binate 
literal l t 


S, 


ISOP - 5, 


AW(S;) + 

N^ISOPS;) 


c 


ad + bda + abde + be 


abc+ acd 


3 


d 


abc + ad+ bda + be 


acd + abde 


3 


d 


abc + acd + ed + abde + be 


ad + bda 


3 



Each row in Table 9 shows 5; and ISOP - S, obtained by Subset{SOP, 
Id for binate literals /, = a,a,b,b,c~c,d<mdd of ISOP. Row 3 (i.e. binate literal b) gives 
the minimum N^) + N BV (SOP - S-J and so is selected. The selected S, is still 
binate, so the second iteration of the outer loop is performed with the 5, assigned to 
SOP; see Table 10. Each row gives Subset(SOP, /,) for binate literals l t = d and d 
in SOP. 



TABLE 10 





SOP=S besr 


= ad+ bda + ed + be 




Binate 
literal 


s t 


ISOP - s t 


N BV iISOP-S t ) 


d 


ad + bda+ be 


ed + abc + acd + abde 


1 


d 


ed+bc 


ad + bda + abc 
+ acd + abde 


3 



The first row (i.e. literal d) of Table 10 gives the lower AW(Si) + 
N BV (SOP - 3), and is selected. The selected $ now is unate and so the process is 
done. We finally obtain g u = ad+bda+bc and 

g b = ed+ abc + acd+ abde . Since Find-Unate-Cube-Set{ISOP) aims 
to reduce both N BV (SJ and N BV (SOP - it tends to make g b highly unate as well. 
Observe that the g b produced in this example is unate for all but one variable (a). 

The Subset(SOP, /,) operation conducted using an M-node ZSBDD 
has a complexity of 0(M). Find-Unate-Cube-Set{ISOP) for an ISOP with TV binate 
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variables repeats the inner loop iV 2 times. Hence the worst case complexity of Find- 
Unate-Cube-Set(ISOP) is OQPM). Compare this with the complexity 0(2 m ) of an 
exact method discussed above; m is usually significantly greater than N and M. 
Thus the presented AND-OR decomposition process can generate highly unate gj 
5 and g b j quite efficiently. 

While embodiments of the invention have been illustrated and 
described, it is not intended that these embodiments illustrate and describe all 
possible forms of the invention. Rather, the words used in the specification are 
words of description rather than limitation, and it is understood that various changes 
10 may be made without departing from the spirit and scope of the invention. 
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